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Abstract 

Quark mass effects are analyzed at high Q 2 in the current fragmentation region of DIS. 
It is found that the linear combination F2 — 2.75F 2 C scales at large Q 2 and small x. We 
obtained a lower bound for the ratio F2/F2 which lies very close to the data from HERA. 
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1 Introduction 



The data on open-charm production in Deep Inelastic Scattering (DIS) from HERA col- 
lider [H, U show that its contribution, F 2 , to the structure function runs up to 40% at 
measured values of x and Q 2 and increases faster than F 2 when x becomes smaller. The 
contribution of b-quarks, F 2 , to the total structure function is of the order 2-3%, as one 
can see from recent measurements of open-beauty production |j| . 

One often believes that, with increasing energy of colliding particles, W, and momen- 
tum transfer squared Q 2 , mass effects become insignificant. However, arguments were 
given in [|], |5| that the difference between the DIS structure functions with open heavy 
quark production in the current fragmentation region and structure functions of the pro- 
cess without heavy flavours scales for large Q 2 , i.e. it depends on the Bjorken variable x 
and heavy quark mass tuq only. This result gives a possibility to obtain model indepen- 
dent (i.e. independent on the concrete choice of the gluon distribution inside the nucleon) 
lower bound for F 2 , which is in an agreement with the experimental data on F 2 [f|, [|. 

In the present work we analyse the problem of mass scale influence on the behaviour 
of physical quantities in the framework of the operator product expansion ||. 

It is convenient to introduce the quantity 

F 2 QQ = F 2 Q /e 2 Q (1) 

instead of F 2 {Q = c,b). Analogously for the light quarks q = u,d,s (that we treat 
massless) we define 

Ff = fj/ej , (2) 

where eQ( g ) are quark electric charges. 

The operator product expansion gives the following expression for the F 2 ®\ 

±F 2 QQ (x,Q 2 ,m 2 Q ) = Cg (^,^f,a s ^ 2 ))®f a ^ 2 ){x} + (3) 



x ^ I* I 1 

,2 

,2\\ ~ e I ..2 



! (\^^,a s (ir)) f Q (tr)[.r} + 



.2\\ ^ f f ,.2\ 



r\2 2 

: C q {^^,a s {^))®f q {^)[x\ 



In Eq. (^|) the quantities Ci are coefficient functions, fi are matrix elements of correspond- 
ing composite operators, which can be identified with distributions of quarks and gluons 
inside the nucleon, and /i is a renormalization scale of the composite operators. Symbol 
(g) means convolution in the variable x, 

a ® b[x] = J J<v)b (j) • (4) 

In Eq. (|3|) we have neglected higher twists contributions. To simplify notations, we do 
not show below a s -dependence of the coefficient functions. 
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It is well known that C{ and fi separately depend on a renormalization scheme. Usu- 
ally, the MS-scheme is more preferable than the MOM-scheme, because of more compli- 
cated calculations in the framework of the latter. 

Nevertheless, the MOM-scheme has some advantages. One of them is the universality 
of the calculational algorithm for coefficient functions in any order of perturbation theory. 
Advantages of use of the scheme with momentum subtraction in the case of heavy quarks 
were noted in 0, where authors proposed a mixed, so called, CWZ renormalization 
scheme. 

The purpose of the present work is to show that effects related to the heavy quark 
mass ttiq, "survive" in the large Q 2 limit. There is a linear combination of the structure 
functions F 2 and F£ that has scaling properties. It appears that this scaling takes place in 
different renormalization schemes. The lower bound for the ratio F^jF^ as the function of 
x at fixed values of Q 2 is also calculated. The results are compared with the experimental 
data. 



2 Asymptotic relations between structure functions 

We are interested in the behaviour of F 2 at large Q 2 and small x. We suppose that the 
production of heavy quarks in this region results from gluons, and heavy quarks are not 
included in the evolution of light quarks and gluons. 
So, we can write 

-F 2 QQ (x,Q 2 ,m 2 Q ) = C g {%^)®f g ^ 2 )[x}+ (5) 

Taking /z 2 = /z 2 ,, where Aqq-q <C /i 2 , ^ Q 2 > an d neglecting "intrinsic charm" ("beauty") 

inside the nucleon, i.e. assuming, that there is a scale /i 2 ,, at which heavy quark distribu- 
tion function is small as compared to the gluon one, we simplify Eq. (||) to the relation 

-F^(x, Q 2 , m 2 Q ) = C g (%, % ® fM) [x] . (6) 

x /i /i 

Let us now define the quantity 

AF 2 = Ff(x, Q 2 ) - F? Q (x, Q 2 , m%) . (7) 

Assuming that at small x the main contribution to the DIS structure function without 
heavy quark production is determined by gluons, i.e. by the formula analogous to (^), we 
find 

V = AC 9 ®/ 9 (^)[x], (8) 

where 

^C 9 = C g (y,^,0)-C g (y,^,^). (9) 
po ro ^0 
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Calculations of the gluon coefficient function in the order 0(a s ) in the MOM-scheme 
were given in M. Using the expression obtained there, we get from (0): 



AC g ^ 



00 



•57T 



(y 2 + (l-y) 2 ) In 



1 + 



(10) 



m%(l + 2y(l-y)) ) 

m Q + ^o2/( 1 -y) J ' 

Thus, the quantity AF 2 tends to the finite limit AF 2 (x, tuq) at large Q 2 . It was shown 
in ||, that this result is not the artifact of the MOM-scheme, although the very expression 
for the quantity A g depends, of course, on the renormalization scheme. 

Starting from the expression for AC g , Eq. (0), it is easy to see, that for y < 0.1 
(small y region is most important in the analysis of the behavior of the structure function 
at x < 1) 

ACg > . (11) 

It follows from the explicit form of C g (y, %) in the MOM-scheme for the massless case 
in the order 0(a s ) for any Q 2 and fi 2 that || 



rn 



(12) 



Then from Eqs. 



we conclude: 



> AFo(x,mn) > . 



(13) 



Inequalities (0) are used below to obtain the lower bound for the ratio F 2 C /F 2 . 



3 Charm contribution to the structure function 

It has been found in the previous section (see Eqs. flS|),(|T0"D), that the difference between 
contributions of light and heavy flavours to the DIS structure function scales for Q 2 — ► oo. 
Taking this result as a basis, it is easy to see, that the following linear combination || 

Z a (x,Q 2 ) = F 2 (x,Q 2 ) + aF^(x,Q 2 ,m 2 c )- (14) 
- (4a + ll)F 2 b (x,Q 2 ,m 2 ) 

has scaling behaviour, with a being an arbitrary constant. 

To cancel the contribution of b-quarks from (14), we have taken a = —2.75. Then we 
obtain the prediction, that the linear combination 

S = F 2 - 2.75F 2 C (15) 

must tend to a function, that depends only on the Bjorken variable x (and the heavy 
quark mass) for Q 2 — > oo. 
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Using the expression for AC g in the leading order in a s , we find, that the above 
difference ( [To]) tends to the scaling limit in the region tuq <C Q 2 in the following way 



x 



^0 



A^(^) 
Ho 



(16) 



+ 



Q 2 



where 



h(y) = -y(l-y)[(2-3yf + 3y 2 } . 



Since m 2 — 7 m 2 > 0, we conclude, that the correction in the expression for E(x, Q 2 ), QI5D, 
is positive and decreasing on Q 2 . 

To compare our results with the recent experimental data, we parametrize _F 2 C according 
to the expression ([IB]) for S(x, Q 2 ,rriQ) and fit the data from HERA collider [||. 

For F 2 we use the parametrization of HI collaboration Q: 



with parameters 



ax b + cx d (l + e^fx) ^ 



lnQ 2 + /ln 2 Q 2 + A 



a 


b 


c 


d 


3.1 


0.76 


0.124 


-0.188 


e 


f 


g 


h 


-2.91 


-0.043 


3.69 


1.4 GeV 2 



For F£ we choose the expression qualitatively coherent to the asymptotic behaviour 
of the quantity £ = F2 — 2.75F£ in the variable Q 2 (|16|) 



F 2 c (x,Q 2 ) 



2.75 



F 2 (x,Q' 



ax b (l 



x) 



l + x 5 Ai n Q2 

Q 



where F2 is defined above (p!8|). Fitting the data from HERA 
130 GeV 2 gives the following values of parameters 



(19) 

for 6.5 GeV 2 < Q 2 < 



a 


b 


c 


9 


h 


0.28 


0.15 


-0.08 


5.00 


1.86 GeV 2 



and x 2 /n.d.f. = 34.6/36 = 0.96. 

In Fig. 1 we show the dependence of the quantity S as a function of Q 2 for two values 
of x, for which we have the array of experimental data, obtained at different Q 2 and x, 
closest to the given values x = 0.01 and x = 0.001. As we see, experimental data are in 
a good agreement with our result on the approach to scaling behaviour from above but 
the very existence of the scaling may be checked at higher Q 2 . 

As has been found in 0, inequalities ( |i~3"D allow us to obtain the following estimation 
for the ratio of structure functions: 



F 2 (x,Q 2 ) 



> 0.4 1 



F 2 (x, m 2 ) 
F 2 (x,Q 2 ) 



(20) 
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It is important to stress, that we did not use any parametrization for F% to obtain the 
inequality, and also, that it does not depend on the behaviour of the gluon distribution. 

Curves represented in Figs. 2, 3 are calculated by using the formula ( ^0|) for two values 
of c-quark mass and are compared with the data of ZEUS collaboration J|. Despite the 
fact that these curves are only lower bounds for the ratio F2/F2, they lay very close to 
the experimental points. 

Our estimations show that lower-bound curves are also in a good agreement with new 
preliminary data of ZEUS collaboration M, including the data for maximal measured 
value Q 2 = 565 GeV 2 . 

4 Conclusions 

In the present work we analyzed quark mass effects in DIS with the help of OPE. It is 
shown by calculating in the leading order that the new scaling takes place in DIS; certain 
linear combination of the DIS structure function and the DIS structure function with 
open charm production scales in the limit of large Q 2 . 

It is found that this specific scaling is in a good agreement with the Q 2 -trend of recent 
experimental data on F| and F 2 obtained at HERA. 

Certainly, higher orders could change the theoretial conclusion about such a scaling. 
At present the first order result can serve as an indication of the existence of an interesting 
physical phenomenon. Influence of higher orders is a subject of our further work. 

We also calculated the lower bound for the ratio -F|/F 2 as a function of x at fixed Q 2 , 
independent on the shape of the gluon distribution in the nucleon, and compared it with 
the data from ZEUS collaboration. This lower bound appears to be quite close to the 
data. 

We are grateful to L.K.Gladilin for providing us with recent results of ZEUS collabo- 
ration on the ratio F 2 C /F 2 . 
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Figure captions 



Fig. 1: Q 2 -dependence of the difference of structure functions at two fixed values of x. 
Solid curves are obtained by fitting the data on F 2 C from [|. The experimental 
points have x closest to the given values x = 0.01 and x = 0.001. 

Fig. 2: The ratio F2/F2 as the function of the variable x at fixed values of Q 2 . Dashed 
curves represent the lower bound for F2/F2 for the c-quark mass m c = 1.7 GeV, 
dotted curves correspond to m c = 1.3 GeV. Experimental points are taken from 0. 

Fig. 3: The same as in Fig. 2, but for other values of Q 2 . 
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Figure 1: 
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Figure 2: 
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Figure 3: 
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